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The physics of magnetic state change in single domain magnetic grains (called Stoner particles) is 
interesting from the fundamental as well as the applied points of view. A change in magnetization 
can be finely tuned with a specific time variation of an externally applied magnetic field. It may also 
occur naturally (without application of a field) at very low temperature with quantum tunneling 
and at higher temperature with thermal excitation. The optimal (usually shortest) time altering 
the magnetisation along with the smallest applied magnetic field are sought in technological appli- 
cations such as high-density reading or writing of information, spintronics, quantum information 
and quantum communication systems. 

This work reviews the magnetization change with a time dependent field and temperature and dis- 
cusses the time it takes to alter the magnetization as a function of the control parameter chosen, 
temperature and material parameters. 

PACS numbers: 51.60.+a, 74.25.Ha, 75.00.00, 75.60.Ej, 75.60.Jk, 75.75.+a 

Keywords: Magnetic properties. Magnetic materials. Hysteresis. Magnetization reversal mechanisms. Mag- 
netic properties of nanostructures 



I. INTRODUCTION 

The effect of time dependent fields on magnetization 
state is important for the reading and writing of infor- 
mation and the monitoring of the magnetization in a 
magnetic material. In the case of magnetic recording, 
when density is increased, the grain size making the 
recording media decreases. If it is small enough, its 
magnetization becomes extremely sensitive to thermal 
energy; it can flip (in the vertical case) or reverse (in 
the horizontal case) by the simple effect of (even small) 
finite temperature perturbation (Brownian fluctuations). 
This adverse effect is called super-paramagnetism that 
traditionally limits longitudinal recording of hard disks 
to densities on the order of 100 Gbits/in 2 . Longitudinal 
refers to the fact the rotation velocity of the disk is par- 
allel to the magnetization orientation. In perpendicular 
recording this limit is ten times higher around several 
Tbits/in 2 . 

Progress in size reduction toward the nanometer paves 
the way to new opportunities in the emerging field of 
spintronics On that scale, we have a wide panel of 
physical effects (e.g. new types of quantum exchange 
between nanometer thick magnetic layers) and the spin 
diffusion length becomes long enough to maintain useful 
spin orientation. Novel nanometric magnetic devices are 
good candidates for use as building blocks of spintronics 
(spin diode or spin transistor) or quantum information 
systems. The latter span quantum information storage 
(Q-bits), quantum computing (quantum logic operations 
like the square root of the NOT operation (y/NOT) or 
the controlled NOT operation (CNOT), the Quantum 
Fast Fourier Transform...), quantum communication 
systems (an example is entanglement which means that 
measurement performed on one system seems to be 



instantaneously influencing other systems related to it) 
or Quantum Metrology. 

It is important to be able to tell how one might be 
able to alter some state with a magnetic field, or how 
it might be affected by temperature as in the recording 
case. It is important to point out that interesting 
quantum phenomena might occur at low temperatures 
(T ~ OK) such as magnetic quantum tunneling; while 
this is beyond the scope of this paper, the reader might 
consult a review such as ref . 0] • 

The time it takes for an effect to take place is also 
important. In recording, given a fixed rotational velocity 
of the hard disk (typically 7200 rpm) the decrease of bit 
length, imposes a faster (higher frequency) process of 
reading (sensing the magnetization orientation) /writing 
(changing the magnetization orientation) of the bit. 
While the shortest time altering the magnetization is 
required in reading/ writing applications, the longest 
time is required in (long-term) storage with protection 
of the data against large magnetic fields that might 
corrupt or even erase the stored information. 

In this paper corresponding to the second part of the 
series on the Stoner-Wohlfarth model, we examine the 
effects of time dependent field and finite temperature on 
a single domain Stoner particle [3|. The time it takes 
for the magnetization to change is also studied with 
temperature and material parameters. 

This paper is organized as follows: in section 2 we 
examine the evolution of magnetization state with a time 
dependent field. In section 3 we discuss the effect of 
temperature on magnetization reversal and we conclude 
in section 4 with the possible extensions and perspectives 



2 



of the SW model. 



II. EQUATION OF MOTION FOR THE 
MAGNETIZATION IN THE PRESENCE OF A 
TIME DEPENDENT FIELD 

Magnetization dynamics is governed by the Landau- 
Lifshitz- Gilbert equation. Since magnetization M is akin 
to angular momentum, we have an evolution equation for 
M similar to angular momentum, the Bloch equation of 
motion dM/dt = JqM x H with 70 the gyromagnetic 
ratio and H the external field (see for instance Kittel 

Extending the Bloch equation to a moment M sub- 
jected to an "effective" field H e and a dissipation term 
describing losses and relaxation processes in the material, 
Landau and Lifshitz (L-L) Q assumed that dissipation 
is accounted for by a coefficient A and introduced a dissi- 
pation non-linear term of the form pj^i M x (M x H e ) 
(where \\M\\ is the modulus of M) controlled by the ef- 
fective magnetic field H e : 



dM 
~dt 



= -70 (M x H e ) 



A70 

IIMII 



M x (M x H e ) (f ) 



L-L define, as in Quantum Field Theory, the effective 
field H e from the functional derivative of the total energy 
with respect to magnetization H e = —SE/6M; hence in 
any magnetic problem the total energy should be the 
starting point whether one is dealing with static or dy- 
namical problems. In the simple Stoner particle case, the 
functional derivative reduces to the gradient with respect 
to M, H e = -dE/dM. 

In order to avoid the divergence problem arising in 
the L-L equation for the large dissipation case (A >> I), 
Gilbert modified the L-L dissipation term by introducing 
a damping term of the form wyM. x (tjt). 

The equation of motion of a magnetic moment in 
presence of damping and effective field is given by the 
Landau-Lifshitz- Gilbert (L-L-G) equation: 



dM 
~~dT 



-7(M x H e ) + ajM x ( 



dM, 

~dT' 



(2) 



where M is the magnetization vector, H e the effective 
field, 7 another gyromagnetic ratio and a the damping 
parameter. 

First of all, we retrieve Bloch equation in the simple 
case of zero damping and effective field H e = H the 
externally applied field. In the static case, the L-L-G 
equation reduces to M x H e = meaning the static equi- 
librium condition is either M / / H e or H e = equivalent 
to the extremum (minimum) condition on the energy as 
discussed previously (see first part of this work) . 
The L-L-G equation conserves ||A<f|| = v 7 M ■ M as seen 
by taking the scalar product on both terms of the RHS of 



eq.H One gets ^jf ■ M = meaning that ||M|| = M s = 
constant, where M s is the saturation magnetization. 
The L-L-G equation seems odd from the mathematical 
point of view since one is used (in systems of ordinary 
differential equations or ODE) to see the first derivative 
term in the LHS only. Here it appears on both sides 
and pushes one to think that the system cannot be han- 
dled by standard mathematical integration tools like Eu- 
ler or Runge-Kutta methods. 

In addition, it is misleading to attempt at solving recur- 
sively the L-L-G equation by substituting repeatedly the 
term ^jf in the RHS of the equation. It is straightfor- 
ward to show that the Landau-Lifshitz (L-L) equation 
is mathematically equivalent to the L-L-G equation by 
taking the cross product of the LHS of eq.[2]with M and 
using M norm conservation {^M- ■ M — 0). We obtain: 

M x ^ = - 7 M x (M x H e ) + a 7 M s 2 (^) (3) 
Substituting M x ^jf in eq. d we get: 



dM 
~dT 



or 



r(M x H e 



07 



1 



-M x [M x H, 



(4) 



It is now a matter of interpreting the coefficients ap- 
pearing in the L-L or the L-L-G equations that will make 
them differ in a given situation. If one identifies 70 as 
pr^2 and A70 as then both equations are same but 
if one insists on keeping 70 as the gyromagnetic ratio or 
confusing dissipation (A) and damping (a) then the equa- 
tions will differ since the factors affecting both terms in 
the RHS are numerically different. 

In addition, the L-L dissipation term goes to zero when 
the damping coefficient goes to infinity making the L-L- 
G equation appear more physically appealing than the 
L-L equation. 

The system of ODE eqs.0]is integrable by standard ex- 
plicit methods, such as Euler or Runge-Kutta (see for in- 
stance ref. [6|]) after expressing the components in Carte- 
sian coordinates. The conservation of the norm is very 
useful during integration (specially in explicit integration 
schemes) to test the accuracy and stability of integration. 
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with the definitions: m 
H e /H K and 5 X 
a{M z /M s ). 



M/M s , rh 



a(M x /M s ),5 y 



dm 

dt 



a(M y /M s ),S z 
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Using order-4 Runge-Kutta (RK4) method (see ref. p|) 
with M along the z-axis as an initial condition, we ap- 
ply at t = a time dependent field making 135° with the 
z-axis (see first part of this work). The 3D response of 
the magnetization in time is depicted in fig. [1] and the z- 
component of m is depicted in fig. [2 Undesirable ringing 
effects (oscillations) in the time variation of m are ob- 
served. They are so because they introduce an unwanted 
delay in magnetization reversal. 

In order to eliminate the ringing effect, we move on 
to another reversal mode called precession switching in 
which the field is applied perpendicularly to the initial 
magnetization and whose action is on until the magneti- 
zation is reversed without displaying any ringing effect. 
The reversal path on the unit sphere is called a ballistic 
path (see fig. [3]) emphasizing its optimality. The sensitiv- 
ity of this process stems from the fact, the field must be 
switched off exactly at the time magnetization reverses 
(see figH]). 



III. EFFECT OF TEMPERATURE ON 
MAGNETIZATION DYNAMICS 

A grain at finite temperature is prone to thermal ex- 
citations that might alter its magnetization state. The 
simplest model describing the effect of temperature on a 
grain is inspired from Chemistry and is called the Neel- 
Arrhenius thermal excitation model. 
At very low temperature, switching may occur by tun- 
neling at a given energy through the energy barrier sep- 
arating two magnetization states corresponding to two 
energy minima (see fig. [5j. This is known as Macro- 
scopic Quantum Tunneling of Magnetization that we will 
not describe here but for which there exist many reviews 
(see for instance ref. Q). 

At finite temperature, the empirical Arrhenius model 
is used to describe the kinetics of a thermally activated 
process. This assumes that an energy barrier hinders 
the forward progress of a chemical reaction. The height 
of this energy barrier is a measure of resistance to the 
reaction. Forward progress of the reaction requires the 
supply of an activation energy to surmount this barrier. 
It has the form [7[: 



r = r exp(A£;/A;BT) 



(G) 



where r is the chemical reaction "inverse rate", To is the 
attempt time to traverse the barrier, AE is the barrier 
height, ks is Boltzmann constant and T is absolute tem- 
perature. 

Drawing an analogy from radioactivity, one might 
view switching as a decay process with a typical proba- 
bility of decay A. Starting from an assembly of grains 
No at t = 0, the number of particles that decay in the 
instant [t,t + dt] is dN — —XN(t)dt. Integrating with 
the initial condition No = N(t = 0) we find that the 
particles that are still present (did not decay or switch) 



is given by N(t) = Noexp(~Xt). This analogy holds if 
switching is treated as an irreversible process like decay. 
This means, switching back to the original value is not 
considered as a valid process. This is the case of data 
recording: If the stored value has changed once, it is no 
longer valid and must be rejected. 

Since the average lifetime is given by r = 1/ A we inter- 
pret the inverse rate as the average lifetime with respect 
to switching. This means that the recorded information 
in a magnetic material (tape, hard disk, floppy etc..) 
stays unaltered for a period of time given by t. We 
infer from this analogy that the probability of switch- 
ing is given by exp(— t/r) and therefore the probability 
of retaining the information (not switching) is given by 
the complementary probability: P(t) = 1 — exp(— t/r) 
with the new interpretation of Arrhenius formula r = 
Toexp(A£/fc B T). 

This decay picture of switching can be recast in a two- 
level model since switching means we have a transition 
from a magnetization state (I) to another (2) as depicted 
in fig. 

Considering a number (normalized) of non-interacting 
grains in state (I) as n\ and the number of grains in 
state (2) as n 2 we may write a kinetic equation (Master 
equation) with typical transition times ti,T2 as: 



(7) 



dn\ n\ ri2 

dt T\ T 2 

Assuming total number (normalized) conservation: 
ni + n 2 = 1 , the solution of this equation is given by: 



± [n>i,o 



n + t 2 



Tl + T 2 



|exp(-t/r) (8) 



where n^o is the initial value of n\ i.e. n^o = n\(t — 
0). It is interesting to note that the decay time r = J 1 ^. 
is the geometric average of t\ and r 2 . 

As a result, we obtain a simple classification of the 
possible magnetic states: 

• We have a blocked state when r >> t i.e. n\ = 

"4,0 Vt. 

• We have a super-paramagnetic state in the opposite 
case t << t leading to m,2 



T1+T2 ' 



Physically, t is of the order of the experimental mea- 
surement time and a blocked state means that no change 
to the system is observed during t. On the other hand, 
when the instrinsic time t << t, the magnetization 
change is so frequent that no well defined state is main- 
tained for a long enough time. Thus the system behaves 
like a paramagnetic system that cannot store informa- 
tion (in a stable and reliable way). Hence the origin of 
the "super-paramagnetic" qualifier. 

When a grain switches we have information storage er- 
rors and the bit error rate (BER) is given by the switching 
probability exp(— t/r). 
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In order to appreciate the meaning of BER and there- 
fore average lifetime and barrier height, suppose we im- 
pose a BER of 10~ 12 . This means one bit is wrong in 
a hard disk of 125 GBytes capacity. Identification of 
BER and switching probability exp(— tjr) means that 
t = 10~ 12 x t. According to eq. [6] and with the as- 
sumptions: to ~ 10~ 9 sec and AE/ksT = 68 we get 
t ~ ir x 10 8 sees which means about 10 years of storage 
(1 year ~ n x 10 7 sees). 



A. Thermal average of the hysteresis loop 

Thermal fluctuations induce random orientations of a 
Stoner particle. If the change of orientation is fast with 
respect to our appreciation of the hysteresis loop, then 
we observe an overall mean behaviour stemming from an 
average hysteresis loop. This average hysteresis loop can 
be calculated with several methods. In ensemble aver- 
aging, one considers a single grain in many orientational 
configurations that is making different angles with the 
magnetic field (taken along the z-direction). With time 
averaging, one considers a single grain undergoing differ- 
ent magnetization cycles while the magnetic field is mak- 
ing different angles with the grain axis. Under the Er- 
godic hypothesis (see ref. p|) these averaging techniques 
should yield the same result. Adopting the ensemble av- 
erage, we ought to find for each angle 4> the minimum 
energy angle 9 and every point on the hysteresis loop is 
made from the average over values of <p. We perform the 
averaging in 3D following the original work of Stoner- 
Wohlfarth despite the fact our previous description was 
intentionally limited to 2D. 

Taking the anisotropy axis along the grain long axis 
with polar angle a and azimuthal angle <f> (see fig. [6]) let 
p(a, (j>) denote the PDF (probability density function) of 
the angles a,<j> be uniform over the domains [0, ^] (see 
note [ll() and [0, 2tt]. Hence, the average loop (being the 
projection of the magnetization over the direction of the 
field) is given by: 



m = cos(# + a) obtaining the equation: 



cos(# + a) 



J 27r d(f> J 2 sin(a)da cos(# + a)p(a, 
/ Q 27r dip J Q 2 sin(a)<ia p(a, <j>) 



(9) 

Since the individual PDF are independent, the joint 
PDF: p{a,4>) = p a (a)p ( j,((j)) is decoupled and since both 
PDF are flat, we get: 



w/2 



cos{9 + a) sin(a)da (10) 



The algorithm is now clear: Sweeping over a we find the 
angle 9 that minimises the energy in order to perform the 
integral. In order to optimize the number of arithmetic 
operations, we rather do the following. We transform the 
minimum equation (as done in the first part) in the form: 
sin(#) cos(0) + hs'm(9 + a) = through the replacement: 



l T,l 



-mcos(2a) ± 



(2m 2 - 1) 
2-y/l - m 2 



sin(2a) 



(11) 



with the plus sign for the upper branch and the minus 
sign for the lower branch. Sweeping over values of m 
since \m\ < 1 allows us to find the corresponding values 
of h from which we keep only the minima energy values 
satisfying the equation: cos(20) + hcos(0 + a) > 0. This 
gives us a table that with proper bookkeeping will help 
us find the average loop. The result of the averaging is 
displayed in the fig. [7] and compared in detail (see the 
figure caption) to the SW work. 



B. Langevin dynamics for the L-L-G equations 

At finite temperature, the deterministic L-L-G equa- 
tion is replaced by the stochastic Langevin equation Q 
governing the evolution of M. The effect of temperature 
is contained in a random additional field r\ (stemming 
from thermal white noise) acting on M: 



dM 
~~dT 



-7o(Mx [H e + V ]) 
Xj0 -M x (M x H e ) 



\M 



(12) 



where the additional magnetic field r\ = (j] x ^y,^z) is 
defined by: 

< m >= 0, < Vi{t)Vj(t') >= 2AS lJ 6(t - t') (13) 

where A is the white noise intensity given by [9j A = 
\k B T/~f M s and = x,y,z). 

Let (9, 4>) be the spherical angles of the orientation of 
the moment M . One may view (9, cj)) as a point on the 
surface of the unit sphere. A statistical ensemble of mo- 
ments with different orientations can be represented by 
a distribution of points over the unit sphere W(9, <fi, t) at 
time t. Conservation of probability leads to a continuity 
equation: 



dW 
~dt 



V- J = 



(14) 



similar to electric charge continuity equation. The 
above is in fact a Fokker-Planck (F-P) partial differ- 
ential equation (PDE) as shown below. The (probabil- 
ity) current density definition J — Wv uses the velocity 
v = j^-^jr of the point (6, <f>) on the sphere, whereas W 
plays the role of a charge density. 

Let us specialize to the case of a single angular de- 
gree of freedom and apply standard methods, fioj to 
write the PDE for the conditional probability density 
P = P(x', t\x, Q). The latter expresses the probability 
density of observing x' = 9{ at time t) given the initial 
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state x = ip>{ at time t = 0). We get the following F-P 
equation: 



dP .. .dP 1 D , ,9 2 P 



(15) 



The "mean first passage time" (time, for the Stoner 
particle, to switch) T(x) satisfies an ODE given by (see 
ref. 0): 



A(x) 



dT{x) 
dx 



-i 



(16) 



In our case, x = iM(^) = - ^§|y where 

-£("0) is the energy (per unit volume) of the Stoner par- 
ticle, E(ip) = — M S H cos(V>) + Keff sin 2 (V'), containing 
Zeeman and effective anisotropy terms. H is the exter- 
nally applied field. Besides ^B(ip) = — 7 2 A yield the 
equation for the "mean first passage time" as: 



A 7o dEjiP) 7 2 A 
M s dip tan(V0 



d^ 2 



-1 



(17) 



This second order ODE can be transformed into a first- 
order equation in v(ip) = dT / dip and integrated once with 
the initial condition v(ip = 0) = 0: 



v(x) = 
with: 



sin(x) 



cxp(-/(x)) 



sing 
7o 2 A 



(18) 



A 



470 FA 



(4HV cos x + (3M S cos 2x + /3M S ) (19) 



2VK e 



V is the volume of the particle and (3 = — — • 
Integrating once again to get T(x) and using the defini- 
tion of the thermal transit time t t h as the value T[x = 0), 
we obtain: 



Uh/c = a 



~a(x+b) 2 pi 

dx T, 27 / dye<v + V (20) 

9o) \ 1 ~ X ) Jx 



The angle 6* maximizes the stationary PDF 
exp[— E(6)/kBT] sin(0) and is also given by the condition 
T(9q) = 0. The coefficients a, b, c are given respectively 
by: 



K, 



eff 



b = 



M, 



k B T ' %K e ff ' 7 Ai4T e // 



(21) 



It is interesting to analyze the results at high and low 
temperatures. In the high temperature limit (a ~ 0); 
we get: tth ~ ca ln(2); whereas at low temperature 
{a{l+b) z » 1), we obtain: 



Hh 



a (1 - b 2 ) (l + b) 



exp[a(l + bf] 



(22) 



Identifying the thermal transit time with r we recover 
in that way the Neel-Arrhenius expression: 

tth ~ T = r o ex p(A£'/fcsT), where the prefactor 

T o = I \/J (1-b 2 ) (T+F)' ( 2 ^) 

and the barrier height AZ? = K e ff(l + b) 2 at low tem- 
peratures. 

In the case of arbitrary temperature, the behaviour of 
the transit time versus a (inverse temperature) for var- 
ious field-anisotropy ratios b is obtained numerically as 
depicted in fig. [8] 

One might be tempted to define the thermal switch- 
ing time directly from the behaviour of the probability 
versus time since the F-P equation provide a means to 
obtain that behaviour. Switching is reached when the 
probability P(w/2,t s ) = 0.5. 

In fig. [9j the time dependence of the probability 
P(w/2,t) is displayed for the field-anisotropy ratio b = 
—0.4 (see ref. [12j) and shows a very quick variation above 
some threshold time if one starts initially from all zero 
values of the probability. 

The results are validated by comparison with the an- 
alytical case in fig. [TU1 Even if the steepness of the nu- 
merical results appear to be weaker than the anatytical 
results, the graph provides a strong support for the ap- 
proximate equivalence of both descriptions. 

Analysis of the thermal switching time versus temper- 
ature has numerous technological consequences. Once 
again, the F-P equation provides this kind of informa- 
tion paving the way to the search of the best materi- 
als/conditions that yield the optimal switching time. 
We perform direct time integration of the F-P equations 
to extract the behaviour of the switching time versus 
temperature. As an illustration, using the same field- 
anisotropy ratio 6 = —0.4 as previously, the inverse 
switching time versus temperature is displayed in fig. 1111 



IV. EXTENSIONS AND PERSPECTIVES OF 
THE STONER- WOHLFARTH MODEL 

The SW model is a macrospin description of magnetic 
systems that is extremely rich from the static, dynamic 
and statistical viewpoints. Despite its numerous limi- 
tations (and of the macrospin approach in general) de- 
scribed in the first part of this work, it remains a valid 
starting point for the useful description and basic un- 
derstanding of many (static and dynamic) problems of 
fundamental and applied magnetism. 

The full 3D counterpart of the SW model as done in 
ref. [HI is an important extension and of great interest. 
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We point anew to the fact the averaging of the hystere- 
sis loop done in section 3 related to thermal effects was 
performed in 3D as in SW work for comparison. 

The extension to uniaxial anisotropics of arbitrary or- 
der (a higher anisotropy is of fourth or sixth order like 
in Cobalt ...) or other forms like biaxial, planar, cubic 
(as in solid Ni or Fe) or of several competing types might 
provide a richer behaviour of the loop versus angle. 

The use of arbitrary non-ellipsoidal shape for the 
grain is also challenging given the occurrence of non- 
uniformity of the magnetization. 



The interaction between grains must also be studied 
and gauged with respect to its role in affecting the switch- 
ing of the magnetization. New types of interactions or 
novel types of exchange between grains or with other ob- 
jects might be exploited in spintronic and quantum de- 
vices. 
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FIG. 1: Trajectory of the magnetization tip on the unit sphere 
for a field applied at t=0 in the y-z plane and making an angle 
of 135 °with the z-axis. The damping is a = 0.1, the field is 
in the yOz plane making an angle of 135 °with z-axis. It is 
applied at t = for 9 nanosecs. Its value is 0.5 Hk- 
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FIG. 2: Variation of the magnetization component m z as a 
function of time. The parameters are the same as in fig. Q] 
The straight thick line indicates conservation of ||iVf|| during 
integration. The dotted line is the variation of the applied 
magnetic field with time. The ringing observed due to damp- 
ing is a major cause of delay in reversal. 
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FIG. 3: Ballistic trajectory of the magnetization tip on the 
unit sphere for a field applied at t=0 in the y-z plane and 
making an angle of 90 °with the z-axis. The damping is small: 
a = 0.001, the field is in the yOz plane making an angle of 
90 °with z-axis. It is applied at t = for 0.12 nanosecs. Its 
value is 1.7 Hk- 




FIG. 4: Variation of the magnetization component m z as a 
function of time. The parameters are the same as in fig. [3] 
The straight thick line indicates conservation of ||JW|| during 
integration. The dotted line is the variation of the applied 
magnetic field with time. No delay in magnetization reversal 
is observed due to the small damping suppressing the ringing. 
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FIG. 5: Energy versus angle showing the barrier AE the 
system has to overcome in order to go from state (1) with 
6\ to state (2) with 62- At low temperature the system can 
tunnel from state (1) to state (2). 
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FIG. 6: System of coordinates displaying the anisotropy axis 
in 3D with the applied magnetic field H along the z-axis and 
the magnetization M all in the same vertical plane indicated 
by dashed lines and making the angle <j> with the xOz plane. 
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FIG. 7: The 3D averaged hysteresis loop looks very much like 
the Stoner-Wohlfarth curve except it is less rounded at the 
approximate switching field values. 
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FIG. 8: Normalized thermal switching time t t h/c versus in- 
verse temperature a for various field-anisotropy ratios b = 
—0.3, —0.4, —0.5, —0.6 as we proceed from left to right. 
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FIG. 9: Switching probability at 9 — tt/2 versus time at a 
fixed temperature for a field-anisotropy ratio b — —0.4 (see 
ref. [12] )■ The ten curves corresponding to ten different in- 
verse temperatures uniformly distributed over the interval 
a — [0 — 50] (see fig. [SJ are indistinguishable. Switching time 
is reached when the probability is equal to 1/2. 
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FIG. 10: Comparison between the analytical formula for the 
normalized thermal switching time t t h/c and the numerical 
integration of the F-P equation versus inverse temperature o. 
The analytical formula leads to a very steep variation with a 
in sharp contrast with the numerical case. 
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FIG. 11: Inverse normalized thermal switching time c/t t h ver- 
sus temperature 1/a for a field-anisotropy ratio b — —0.4. 



